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Abstract 

We study the formation, evolution and structure of dissipative fronts pro¬ 
duced by overtaking collisions of relativistic streams, with emphasis on 
strongly magnetized flows. The evolution of the system is followed using 
analytical approach in the simple wave regime and numerical simulations in 
the non-simple wave regime, until a steady-state is reached. The steady-state 
structure of the front is then examined by solving the appropriate jump condi¬ 
tions. The conversion of magnetic energy into kinetic energy is parametrized 
in terms of the Alfven 4-velocity inside the front. The implications for gamma- 
ray jets are briefly discussed. 
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I. INTRODUCTION 

Highly magnetized, relativistic outflows are believed to play an important role in a variety 
of astrophysical systems, e.g., pulsar winds (Kennel & Coroniti 1984; Gallant, et al. 1992), 
Galactic and extragalactic superluminal sources (Burns & Lovelace 1982; Begelman & Li, 
1994; Levinson & Blandford 1996), gamma-ray bursts (Usov 1992; Waxman 1995; Thompson 
& Duncan 1993; Smolsky & Usov 1996), and others. Such flows are likely to be unsteady, 
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giving rise to formation of shocks which can lead ultimately to the dissipation of some fraction 
of the bulk energy, either in the form of thermal energy or, through shock acceleration 
(Blandford & Eichler 1987) or magnetic reconnection, high energy particles with nonthermal 
distributions. As demonstrated by Romanova & Lovelace (1997) recently, the creation of 
dissipative fronts in Poynting flux jets may give rise to 7 -ray flares as often seen in blazars. 

One important difference between unmagnetized and highly magnetized flows is that, in 
the former case the sound speed is limited to c/\/3, whereas in the latter case the Lorentz 
factor associated with the speed of propagation of a disturbance - the magnetosonic speed - 
can be much larger. Consequently, the conditions required for the formation of strong shocks 
are expected to be more stringent in highly magnetized flows. Even when the upstream flow 
is supermagnetosonic, so that strong shocks are created, only a relatively small fraction of 
the incoming flow energy goes to heat the downstream region; the remainder is used up to 
slightly compress the magnetic field downstream (Kennel & Coroniti 1984). The situation 
might be markedly different however in the presence of rapid magnetic field dissipation 
behind the shock. 

Motivated by the above considerations, we study in this paper the formation, evolution, 
and structure of fronts produced by unsteady injection of relativistic streams, with particular 
emphasis on flows in which the energy flux is dominated by electromagnetic fields (commonly 
referred to as high sigma flows). 


II. THEORY 


A. Basic equations 

The stress-energy tensor of a perfectly conducting fluid can be written in terms of u a , 
the fluid 4-velocity, and the magnetic field 4-vector b a = F a/3 up, where F a13 is the dual 
electromagnetic stress-tensor, as 

T af) = rh*u a U P + p*g af) - b a b?. (1) 

Here, g a g denotes the metric tensor, r is the fluid rest mass density, h* = h + b 2 /r and 
p* = p + b 2 / 2 , where p is the proper gas pressure, h the specific enthalpy and b 2 = b a b a . 
Under the assumption of ideal MHD, Ohm’s law reduces to the condition F a ^u^ = 0. The 
MHD equations can then be written in a divergence form as (van Putten, 1993) 

W a T a0 = 0 , 

V a (ru a ) = 0, (2) 

W a {u [a 0 + g a9 c) = 0, 

where c := b a u a = 0 is a conserved constraint. Note that \[A^b a represents now the magnetic 
field measured in the fluid rest frame, the absolute value of which we denote here by B. The 
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above set of equations must be supplemented by an equation of state for the gas. In regions 
of adiabatic flow where u"V cv S' = 0, S being the entropy, the specific enthalpy is given by 
h — 1 + / dp/r, with p = p(r). 

Using eq. (1), the energy flux carried by the flow can be written in the form 

T = r/rfv(l + u\) = rh'y^v, (3) 

where v is the fluid 3-velocity, u \ = b a /Vrh, is the Alfven 4-vclocity with respect to the 
fluid rest frame, and 7 a is the corresponding Lorentz factor. 


B. Simple waves 


To simplify the analysis, we shall consider in the following ID flow with velocity u a = 
( 7 , u, 0, 0), a polytropic equation of state, p(r) = Kr r , and magnetic held perpendicular 
to the fluid 3-velocity. Using Maxwell equations together with the continuity equation it 
is readily shown that B/r is conserved along stream lines (de Hoffmann & Teller 1958). 
Consequently, the proper magnetic pressure can be expressed as, ps = ftr 2 , with k = 
B 2 /8irr 2 . Thus, the problem reduces essentially to that of a ID hydrodynamic how with an 
effective equation of state p*(r ) = Kr T + nr 2 . Next, we define 

2 1 dp* r/lr r_1 + 2 nr 

s h* dr 1 + J_/y r r - 1 -f 2 nr 


Evidently, in the absence of magnetic helds, viz., k — 0, c s is simply the sound speed. In the 
presence of very strong helds, such that the gas pressure is much smaller than the magnetic 
held pressure, c s reduces to the Alfven speed Ua/xa- If the entropy and k are taken to be 
constant throughout the how, one can define the variables da = c s dr/r and A = skill -1 (u). 
Equations (2) then yield (van Putten 1991) 


dx± 

dt 


± A±(a, A) 


dx± 

dx 


0, 


( 5 ) 


with x± — a ± A, and A±(a, A) = ( c s ± v)/(l ± c s v ), where v — u/ 7 is the huid 3-velocity. 
Equation (5) implies that y + and y_ are constant along the characteristics dx/dt = A + and 
dx/dt = A_, respectively. Note that A± is the relativistic sum and difference of c s and v. In 
the case of simple waves one of the Riemann invariants is a constant. Let us suppose that 

= const. Then, A + depends only on the remaining invariant, y + , and it is clear from 
eq. (5) that x+ is constant along straight lines with slopes A + , implying that A + is simply 
the speed of propagation of a disturbance in the observer frame. For weak disturbances for 
which v —> 0, we obtain A + —> c s , verifying that c s is indeed the speed of propagation of the 
disturbance in the fluid frame. 

For illustration, we consider unsteady injection of ID, highly magnetized fluid with the 
following initial and boundary conditions: 
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A (x,t) = X a = const ; 
X(x = 0, t) — A b it) = A 0 + t/r ; 


t < 0, 
t > 0. 


( 6 ) 


The general solution is then given by 

where the subscript (+) has been dropped to simplify the notation. In order to solve for 
t', we must determine A(A). Since we are merely interested in high sigma flows, we can, to 
a good approximation, neglect the gas pressure. This corresponds to the limit K = 0 in 
equation (4). We then obtain A(A) = tanh(3A/2 — A 0 /2 + A Ao ), where sinh(A Ao ) = ua 0 is 
the Alfven 4-velocity of the initial flow (at t < 0). On substituting the later expression for 
A(A) into eq. (7), we obtain a transcendental equation for the retarded time t', which can be 
solved numerically. The solution thereby obtained describes a steepening wave propagating 
at a velocity tanh(A 0 + Aa 0 ), the relativistic sum of va 0 and v. This wave will steepen into a 
shock after time t s at some distance x s from the injection point. The point (t s ,x s ) is a pole 
of the function A (x,t), and can be found from the equation dX/dx = oo. The result is 


t s 


T 

3 


sinh(2A 0 + 2A Ao ); 


x s 


ct s tanh(A 0 + A Ao ). 


( 8 ) 


If we associate the flow acceleration time r with the size of the injection region, R , (this might 
be, for example, the radius of a black hole or the inner boundary of an accretion disk from 
which a jet is ejected), then in the ultra-relativistic limit eq. (8) gives, x s = -R7 2 7 A /3. Note 
that in this limit 7 o7a is just the Lorentz factor associated with the speed of propagation 
of the disturbance with respect to the injection frame, so that ct( 7 g 7 A ) 2 is roughly the 
distance from the injection point at which two disturbances separated by a time r collide. 
This constraint on the location at which the front is created may have important implications 
for the characteristics of the high-energy emission from relativistic jets (Levinson 1996). In 
the case of the Crab pulsar the above criterion (although the different geometry is expected 
to somewhat alter this result) implies that for R ~ 10 6 cm any fluctuations of the pulsar 
wind are not likely to steepen into traveling shocks upstream the standing shock (which 
is produced as a result of the interaction of the pulsar wind with the confining, slowly 
expanding SNR, and is located at a radius of about 0.1 pc), unless 7 q 7a < 10 5 . 


C. The evolution and structure of adiabatic fronts 

The simple wave solution breaks down after a shock is formed (i.e., for t > t s ). To 
study the evolution of the front in the non-simple wave regime, we have solved equations 
(2) numerically. A detailed account of the numerical method used is given in van Putten 
(1993). To test the code, we have compared the numerical results with the analytic solution 
obtained above in the simple wave regime. The agreement is found to be good. 
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An example is shown in fig. 1 . As seen from this example, the collision of the initial and 
injected fluids leads ultimately to the creation of a front consisting of two shocks moving away 
from each other, and a contact discontinuity across which the total pressure is continuous. 
The front propagates at a speed (i.e., the speed of the contact discontinuity) intermediate 
between that of the initial and injected flows, and expands at a rate proportional to the 
relative velocity of the two shocks. Consequently, there is a net energy flow into the front, 
giving rise to a pressure buildup inside. After the injection is terminated, the front reaches 
a steady-state whereby the net energy flow into the front is balanced by adiabatic cooling, 
owing to the expansion of the front. It is then convenient to explore the steady-state structure 
of the front using the appropriate jump conditions. Now, in the ideal MHD limit the specific 
magnetic field, B/r, must be continuous across the shock. This will no longer be true in the 
presence of magnetic field dissipation, which may result from e.g., magnetic reconnection in 
the front, (this would require magnetic field topology more complicated than that invoked 
above). A self consistent treatment of dissipative fronts is beyond the scope of this Letter. 
In order to explore the effect of magnetic field dissipation on the parameters of the shocked 
fluid, we suppose that some fraction of the magnetic field energy is being converted into 
particle energy inside the front by some unspecified mechanism. We can then treat the 
magnetic field behind each shock as a free parameter which lies in the range between zero 
(although the magnetic field is not likely to drop well below the equipartition value) and the 
maximum value obtained in the ideal MHD limit. 

The problem can be most easily solved in the rest frame of the front. In this frame the 
jump conditions decouple and reduce essentially to those of two single, independent shocks. 
Let subscript -(+) refers to fluid quantities leftward (rightward) of the contact discontinuity, 
and denote quantities inside the front by subscript /. The jump conditions at the right (+) 
and left (-) shocks can then be written in the front frame as, 


r/± = 7±(1 “ v±/v a ±)r±, 

r±h± 7±(1 - v±/v s ±) -p± = r f ±h* f± - p* f± , (9) 

r ±h±l±v±{v± - v s ±) + p± = P*f±- 


Here v s ± denote the shocks 3-velocity, v± are the 3-velocities of the fluids outside the front 
(i.e., ahead of the shocks) with respect to the front frame, and r )± are the corresponding 
Lorentz factors. Since we are merely interested in the strong shock limit, we shall restrict 
our analysis, in what follows, to the ultrarelativistic case, viz., » 1. In this limit we can 
approximate the enthalpy inside the front as, hf± ~ [T/±/(T/_)_ — 1 )](pj ± /r/±). Assuming 
that the fluids upstream the shocks are cold (i.e., p± = 0), eqs. (9) yield for the shock 
velocities, 


(2 + u 2 Af± )Yf± - 2 
2 + T f ±u 2 Af± 


( 10 ) 


and for the front pressure and temperature to the right (left) of the contact discontinuity, 
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Pt± 4 r /±(r /± - i)(i + «L)(i + U %±) 

r ±c 2 7± (2 + r /± ^ /± )(2r /± -2 + r /± «^ /± )’ 

kT f± = Pf± _ 2(r /± -i)(i + ^ ± ) 

me 2 7± r/± c 2 7± 2 + r /± ^ /± 

Note that the shocks velocities given by eq. (10) are independent of the magnetic held 
outside the front. 

In the ideal MHD limit, the requirement that the specific magnetic held should be con¬ 
tinuous across the shock implies that the Alfven 4-velocities ahead and behind the shocks 
are related by, 


Vf±u\f ± + [2(r /± - 1) + (r /± - 4 )u 2 A± \u 2 Af± - 2u 2 A ± - 0 , (12) 


where eqs. (10) and (11) have been used in deriving eq. (12). In the weakly (u A ± « 1) 
and strongly (ua± » 1 ) magnetized cases eq. ( 12 ) simplifies to ( ua/±/ua ±) 2 = l/(r/± — 1 ) 
and («a/±/«a ±) 2 = (4 — T/±)/rj-t, respectively. On substituting the latter results into 
eqs. (10) and (11) we recover the shock velocity and downstream temperature of a perfectly 
conducting fluid, obtained earlier by e.g., Gallant et al. (1992); v s ± = ±(r/± — 1), kTf± = 
(r f± — l)mc 2 7 ±, in the weakly magnetized case, and (u s ±/ua ±) 2 = (4 — T f ±)/(8 — 41" f±), 
kTf± = [(2r /± — 2)/(4 — r/-i-)]mc 2 7 ±, in the highly magnetized case. We see that in the 
strong shock limit considered here, a perfectly conducting front is heated to relativistic 
temperatures for any value of the magnetic held. However, while in the case of weakly 
magnetized streams all of the bulk energy per particle incident into the front as measured in 
the front frame, 7 ±mc 2 , is converted into random energy (the average energy per particle in 
a relativistic gas is kT/\T /± — 1]), in the highly magnetized case only a small fraction of the 
total bulk energy carried by the how outside the front, specifically the fraction carried by 
the particles ~ 7 ( see ec f [3]), is thermalized; the remainder of the energy is used up to 
compress the magnetic held inside the front. Consequently, collisions of highly magnetized 
hows in the ideal MHD limit lead to the formation of magnetically dominated fronts. 

The situation could be markedly different if the front is highly dissipative. Then, u a /± 
can deviate substantially from the maximum value given by eq. ( 12 ). As seen from eq. (11), 
the energy per particle, kTf±/ (T f± — 1), increases as uaj± decreases, ultimately approaching 
m( ?l±l\± in the limit of rapid magnetic held dissipation, i.e., Ua f± < 1. (We have assumed 
that the front is adiabatic. It is conceivable that it will become radiative at temperatures 
lower than that derived above if the synchrotron or inverse Compton cooling time is suf¬ 
ficiently short. We shall not consider this possibility in this Letter.) The shock velocities 
approach v s ± = ±(T f± — 1 ) in this limit. 

Now, the velocity of the contact discontinuity and, hence, the velocities upstream of each 
shock with respect to the front frame, u±, are unknown a priori. The requirement that the 
total pressure inside the front is continuous across the contact discontinuity surface yields, 
in the strong shock limit, the relation, 


r„T/-(:i + u\_){ 1 + u\ f _){ 2 + T f+ u 2 Af+ ) 
r + T /+( 1 + Ua +)( 1 + m a /+)(2 + r 
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Note that in the case of a symmetric front, which prevails when r_ = r+, = p + , and 

ua- = ua+, eq. (13) gives r/ = 1 as it should of course. 

To complete the solution, we transform to the frame in which the fluid is injected (here¬ 
after referred to as injection frame; this could be e.g., the rest frame of a star or a black hole 
from which a jet is emerged), and in which the velocities of the initial and injected flows, 
denoted here by u' + and u'_, are specified. The proper mass density, pressure, and magnetic 
field are of course relativistic invariants and have the same values in the injection frame as in 
the front frame . We consider the case u'_ » u' + > 0, namely both the initial and injected 
fluids propagate in the same direction. On transforming to the injection frame we obtain, 

7± = 7±7c - u±u c , (14) 

with u c being the front 4-velocity with respect to the injection frame. Equations (13) and 
(14) can be solved now to yield y c and y_. One finds 



For initial fluid at rest (u' + = 0) the above solution simplifies to y_ ~ (y(_/2y) 1 / 2 , and for 
relativistic initial fluid to y_ ~ (qL/dr^y^) 1 / 2 . On substituting the latter results into eqs. 
(10) and (11) we obtain an expression for the front quantities in terms of u'±. Evidently, for 
roughly symmetric fronts (77 order unity) the front temperature scales roughly as ( 7 k/Gy^) 1 / 2 . 
In the asymmetric case, one side of the front will be heated to much higher temperatures 
than the other. Heat conduction across the contact discontinuity surface may then become 
important if the magnetic field is not completely transverse there. 

III. DISCUSSION 

A highly conducting front produced by overtaking collisions of strongly magnetized, 
relativistic flows should be magnetically dominated, and is not an efficient converter of bulk 
energy into random energy (i.e., thermal energy or nonthermal distributions of accelerated 
particles). Much higher temperatures can be attained in the presence of rapid magnetic field 
dissipation inside the front, that may arise from magnetic reconnection or plasma instabilities 
near the shocks. A fraction of the dissipated field energy is likely to be radiated away in the 
form of electromagnetic emission with relatively high efficiency, owing to synchrotron and 
inverse Compton cooling of the energetic particles. The latter process can be very effective 
in the presence of background radiation with roughly isotropic distribution, as is believed 
to be the case for Galactic and extragalactic jets (Dernier & Schlickheiser 1993; Sikora et 
al. 1994; Ghisellini & Madau 1996; Levinson & Blandford 1996). The radiative friction 
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may also provide a mean for dissipating the magnetic field. Such fronts may explain the 
transient emission observed in many of these sources (Romanova & Lovelace 1997). However, 
in order to be efficient radiator, the front must be created at not too large a distance from 
the putative black hole, depending on the intensity and spectrum of the external radiation, 
thereby imposing a constraint on the magnetosonic speed (see discussion following eq. 8). 
In the case of streams consisting of electron-positron plasma, the front may become loaded 
with e ± pairs as a result of effective pair production. Frequent creation of such fronts may 
lead to a transition from magnetic to particle dominated jet (Levinson 1996), as seems to be 
suggested by observations (Blandford & Levinson 1995). Again, the location at which the 
front is created might be important; the fronts must form above the cooling radius, which 
depends upon the intensity and spectrum of the background radiation, in order to avoid 
catastrophic X-ray production and the disruption of the jet. 

We thank R.D. Blandford for useful comments. A.L acknowledges support from NSF 
grants AST 91-19475 and AST 93-15375, NASA grant NAG 5-3097, and Alon Fellowship. 
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Figure Captions 

Figure 1. Three epochs are shown in two examples of a ID simulation of the formation of 
a relativistically hot, magnetized front in which the magnetic field is everywhere transverse 
to the fluid flow. Displayed are the temperature, the specific magnetic field, and the fluid 
3-velocity in the reference frame co-moving with the initial flow (ahead of the front). A 
burst of new outflow from the left is modeled by the initial 4-velocity distribution u x (x) = 
sinh{l + tanh(10^ 2 ^)}, where x is the coordinate along the outflow, x Q represents the 
injection point and L is the length of the displayed computational grid. The initial density 
(r = 1), pressure (P = 0.025) and Alfven 4-velocity [a a = 0.5 in (a) and ua = 1 in (b)] are 
taken to be uniform. As seen, the higher Alfven velocity in (b) results in a broader front. 
Note the gradual steepening of the leading shock front, which is more delayed in (b). The 
left shock becomes steady at a very early stage in both examples shown. 
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